Magnetic properties of interacting, disordered electron systems in d = 2 dimensions 
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We compute the magnetic susceptibilities of interacting electrons in the presence of disorder on 
a two-dimensional square lattice by means of quantum Monte Carlo simulations. Clear evidence is 
found that at sufficiently low temperatures disorder can lead to an enhancement of the ferromagnetic 
susceptibility. We show that it is not related to the transition from a metal to an Anderson insulator 
in two dimensions, but is a rather general low temperature property of interacting, disordered 
electronic systems. 



I. INTRODUCTION 

The interplay of disorder and Coulomb repulsion be- 
tween electrons raises many fundamental questions not 
only in solid state physics^— but also in the physics of 
cold atoms£~— Electron repulsion can lead to a Mott 
insulator^ a description of which must invoke many- 
body physics. On the other hand, presence of disor- 
der in a system of non-interacting electrons can cause 
Anderson localization^^ The simultaneous presence of 
disorder and electronic interaction is known to give rise 
to new, subtle many-body phenomenal— An impor- 
tant question in the physics of interacting and disor- 
dered electrons is the stability of magnetic order. The 
effect of disorder has on magnetic long-range order has 
been explored by various analytica l 10 ' 11 and numerical 
methods^ - — A microscopic lattice model to study these 
effects is the Anderson-Hubbard model, which has been 
explored using quantum Monte Carlo (QMC)j^ 2 - dynam- 
ical mean-field theory (DMFT)^^ and Hartree-Fock 
calculations! 16 ' 17 One of the questions that have received 
particular attention is whether the long-range antiferro- 
magnctic order (AFLRO), that occurs at half- filling in 
the clean limit, is stable when disorder is present. In- 
creasing disorder has gencrically been found to be detri- 
mental for AFLRO 

Magnetic instabilities and the formation of local mo- 
ments have also been investigated in high mobility semi- 
conductor heterostructuresJ^ Some experiments have 
found the ferromagnetic (uniform) susceptibility to be- 
have critically near the MIT^ while other experiments 
only found a strong enhancement i 20 ' 2 1 

Analytical work has predicted critical behavior of the 
ferromagnetic susceptibility as the metal-insulator tran- 
sition is approached from either side of the transition, 
and the physics has been attributed to either the forma- 
tion of local moment a 22 ' 23 or to the vanishing of the spin 
diffusion constant cither at or near the (charge)transport 
critical point i 10 ' 24 

QMC investigations in two dimensions have also sug- 
gested the existence of a quantum critical point (QCP) 
between a metallic and an (Anderson) insulating phase 



when the disorder strength is variedi^— Magnetic prop- 
erties both near and away from the QCP have received 
less attention. Denteneer et aL— found an increase in 
the ferromagnetic susceptibility both sides of the metal- 
insulator transition and this behaviour was attributed to 
the formation of local moments ! 22 ' 23 

In this paper, we report results on the magnetic prop- 
erties of the Anderson-Hubbard model in two dimensions 
obtained through extensive QMC simulations. We con- 
centrate on the ferromagnetic (FM) and the antifcrro- 
magnctic (AF) susceptibilities. There are two main re- 
sults: (i) We find that, at sufficiently low temperatures, 
an increase in the disorder increases the FM susceptibil- 
ity, and (ii) the enhancement of the FM susceptibility is 
related neither to the MIT, nor can it be explained by 
the formation of local moments. In contrast, increasing 
disorder always reduces the AF susceptibility. 



II. CORRELATION FUNCTIONS FOR THE 
ANDERSON-HUBBARD MODEL IN d = 2 

Our investigation of interacting electrons in the pres- 
ence of disorder is based on the Anderson-Hubbard 
Hamiltonian on a square lattice 



H = T{e,J 



Here 



(1) 



(2) 



<ij >a 



is the single-electron part where c ia (c.; CT ) are fermion 
creation (annihilation) operators for site Ri and spin cr, 

Uic = c\ a Cia is the operator for the local density, /i de- 
notes the chemical potential, and t is the hopping am- 
plitude for electrons between nearest neighbor sites. The 
local energies q are random variables which are sam- 
pled uniformly from the interval [—A/2, A/2]; hence the 
width A characterizes the strength of the disorder. The 
interaction is taken to be repulsive (U > 0). 
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The model is solved numerically using detcrminantal 
QMC (DQMC)3 The hopping integral t sets the unit of 
energy (we set S=fcs=l) and the simulation now contains 
three independent energy-scales: the disorder strength 
A, the interaction strength U, and the temperature T. 
The average filling n is tuned by the chemical potential 
\i. 

To extract the magnetic behaviour, we focus on on the 
real-space equal-time spin-spin correlation function^ de- 
fined as 



C(r) = (S 2 (R,+r)S 2 (R 4 )>, 



(3) 



where S z (Ri) = — mi- C(r) measures the extent to 
which the z-component of a spin on site is correlated 
to the z-component of another spin at a distance r. Al- 
though we have used the operator S z in Eq. [3] to define 
C(r), we explicitly confirmed that our calculations have 
the full SU(2) invariance of the problem. 

From C(r), we define the spin structure factor 5(q) as 
the Fourier transform: 



S(q)=$>*"C(r). 



(4) 



The q-dependent susceptibility is obtained from x(q) — 
/3S{q)^- The wave-vectors q = (0,0) yield the FM sus- 
ceptibility \f, and q = (w, ir) yield the AF susceptibility 
Xaf, respectively. The lattice spacing a is assumed to be 
unity. All plots shown are for a 10 x 10 square lattice. 
We have confirmed that the conclusions are unchanged 
for computations on 8 x 8 and 12 x 12 square lattices by 
explicit simulations. 

There are two sources of statistical error in this analy- 
sis: one is due to the QMC simulations, the other arises 
from disorder averaging. For all parameter sets stud- 
ied here, the intrinsic QMC error for any given disorder 
realization is much smaller than the error arising from 
different disorder realizations. 



III. RESULTS 

We start with the non- interacting system. In Figs.QJa) 
and (b) we show the dependence of the inverse AF and 
inverse FM susceptibilities, respectively, on T for U = 
and various disorder strengths. Without interaction, an 
increase in disorder always reduces both the AF and FM 
susceptibilities. Namely, increasing disorder creates more 
fluctuations in the site-energy landscape and, since U = 
0, two electrons of opposite spins will tend to occupy the 
lowest energy sites, thereby reducing the effective local 
moment and destroying magnetism. The inset plots the 
square of the local moment as a function of temperature, 
defined as (m 2 z i ) = C(r = 0) = (n*) — 2(71^^44.). 

Now we include the interaction U. In Fig. [21 we plot 
the inverse AF susceptibility for two different (U, n) pa- 
rameter sets. In all cases, 1 /xaf decreases linearly with 
decreasing T. It is clearly seen that the plot for ev- 
ery value of A can be extrapolated to zero at T = 0, 
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FIG. 1: (Color online) (a) Inverse AF susceptibility and, 
(b) inverse ferromagnetic susceptibility computed for a 10 
x 10 square lattice at quarter-filling (n = 0.5) for U = 0. 
The susceptibility data is averaged over 10 disorder realiza- 
tions for the highest temperatures (T = 1.0,0.5,0.333,0.25), 
80 disorder realizations for intermediate temperatures (T = 
0.2,0.167), 100 disorder realizations for low temperatures 
(T — 0.125,0.1) and 120 disorder realizations for the lowest 
temperature (T = 0.0833). Inset: Square of the local mag- 
netic moment vs T from the same datasets. The horizontal 
line is the local moment squared (= 0.375) when A = 0. 



which signifies that there will be no AF ordering. This 
behaviour is expected for the Anderson-Hubbard model 
far away from half-filling (n = 1). Furthermore, the 
AF susceptibility decreases monotonically as the disor- 
der strength is increased. 

This is not surprising since AFLRO will be hindered by 
randomness. Namely, electrons inevitably find it difficult 
to hop (and thus gain kinetic energy) when the energies 
on two neighbouring sites are very different. Ferromag- 
netism, on the other hand, stems from the need to avoid 
strong Coulomb repulsion locally. In Fig.[3l we plot 1/xf 
vs. T for the same two parameter sets as in Fig. At 
high temperatures Xf decreases with increasing disorder 
strength. However, at low temperatures, the opposite 
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FIG. 2: (Color online) Inverse of the AF susceptibility for the 
parameter sets (a) U = 2, n = 0.3 and (b) U = 3, n = 0.5. 
The number of disorder realizations is the same as in Fig. 
Q] Similar behaviour is observed for all other parameter sets 
investigated. 
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FIG. 3: (Color online) Inverse of the ferromagnetic suscep- 
tibility for the parameter sets (a) U = 2, n = 0.3 and (b) 
U = 3, n = 0.5. At sufficiently low T (T less than about 0.3 
in (a)) xf increases with increasing disorder strength. 



behaviour is seen to occur. We point out the difference 
between Fig. [3] and the non-interacting case in FigQJb) • 
Thus it is clear that the enhancement of xf with increas- 
ing disorder at low temperatures takes place only in the 
simultaneous presence of interaction and disorder. Also, 
this behaviour is ubiquitous in all parameter sets we have 
investigated. Even at weak coupling (e.g. U = 1), the 
ferromagnetic susceptibility rises with increasing disorder 
at lower temperatures. 

To investigate if the enhancement of the ferromag- 
netic susceptibility is related to the metal-insulator tran- 
sition (MIT) in d — 2, we plot \f in Fig- HI against the 
disorder strength for U = 2, n = 0.3. For these pa- 
rameters, the d = 2 Anderson-Hubbard model shows a 
metal- (Anderson) insulator quantum phase transition at 
A c = 7.8jSS It is clear from Fig. 0]that xf does not show 
any enhancement close to the MIT. This holds for all 
parameter sets we have investigated. 

In Fig. [5]wc plot the average squared moment vs. tem- 
perature for different values of the disorder strength for 
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FIG. 4: (Color online) Ferromagnetic susceptibility as a func- 
tion of disorder strength from high (bottom) to low (top) 
temperatures for U = 2, n = 0.3. 
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FIG. 5: (Color online) Square of the average moment vs. tem- 
perature for different disorder strengths for U = 2, n = 0.3. 
The number of disorder realizations is the same as Fig. [T] The 
magnetic moment continues to monotonically decrease as dis- 
order is increased for all T. The temperature values where we 
observe a sudden change in the slope of {ml) for all disorder 
strengths correspond closely to the temperature values which 
demarcate high and low temperature regimes in Fig. [HJa) . 



in randomness can increase the tendency towards ferro- 
magnetism. This tendency towards ferromagnetism is 
present for all interaction strengths we have investigated, 
although, as expected, stronger interaction increases the 
ferromagnetic susceptibility. We have shown that this en- 
hancement in ferromagnetic suscpetibility is not due to 
(i) the formation of local moments and (ii) the Stoner in- 
stability. It is a ubiquitous low-temperature phenomenon 
in the simultaneous presence of disorder and interaction. 
In particular, this is not related to the phase transi- 
tion between the metallic and Anderson-insulating phase. 
Rather, the enhancement is seen in both the metallic and 
insulating phases. We have also established that at low 
temperatures, the response of the spin sector to an in- 
crease in the disorder strength is markedly different to 
that of the charge sector. Similar effects are also seen in 
calculations using DMFT.— Therefore Anderson local- 
ization cannot explain these effects since it is absent in 
d = oo where DMFT is exact. The details of the DMFT 
calculations will be addressed in a future publication. 



U = 2, n = 0.3. As the disorder strength is increased, the 
magnetic moment is seen to decrease at all temperatures. 
Thus the enhancement of the ferromagnetic susceptibility 
at low temperatures cannot be explained by an increase 
in the number of magnetic moments with disorder. We 
also emphasize that the low temperature increase in xf 
is different from the Stoner criterion for ferromagnetism, 
UN(Ep) > 1, where N(Ep) is the electron density of 
states at the Fermi level. Since increasing disorder (at 
fixed U) inevitably leads to loss of spectral density at the 
Fermi level, the Stoner criterion predicts that ferromag- 
netic tendencies are acted against by increasing disorder. 

In Fig. [SI we plot the charge susceptibility \c = 
d(n)/dfi vs. T for different disorder strengths for [7 = 2, 
n = 0.3. The charge susceptibility decreases monoton- 
ically with increasing disorder strength for all temper- 
atures investigated. Since the charge susceptibility is 
expected to be proportional to the density of states at 
the Fermi level (at fixed interaction strength and elec- 
tron density), we infer that N(Ep) is decreasing with 
disorder. The charge susceptibility does not show any 
discernible change in behaviour with disorder at high or 
low temperatures, as in the case of the dc conductivity 
Cdc^ but in contrast with xf and (m^), thus implying 
that at low temperatures the charge and spin sectors dis- 
play different behaviour in the simultaneous presence of 
disorder and interaction. 



IV. DISCUSSION 

In conclusion, our results show that in a strongly cor- 
related system at low enough temperatures, an increase 




T 



FIG. 6: (Color online) The charge susceptibility \c vs. T for 
different disorder strengths for U = 2, n = 0.3. The num- 
ber of disorder realizations is the same as Fig. Q] Note that 
Xc is related to the compressibility k through k — Xc/(n} 2 . 
Since (n) is held constant for all temperatures and disorder 
strengths, k and Xc are proportional. 
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